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Abstract 

We study the phase transition in the abelian lattice gauge theory using the Wilson-Polyakov line as the 
order parameter. The Wilson-Polyakov line remains very small at strong coupling and becomes non-zero 
at weak coupling, signalling a confinement to deconfinement phase transition. The decondensation of 
monopole loops is responsible for this phase transition. A finite size scaling analysis of the susceptibility 
of the Wilson line gives a ratio for j/u which is quite close to the corresponding value in the 3-d planar 
model. A scaling behaviour for the monopole loop distribution function is also established at the point of 
the second order phase transition. A measurement of the plaquette susceptibility at the transition point 
shows that it does not scale with the four dimensional volume as is expected of a first order bulk transition. 
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I. INTRODUCTION 



Lattice gauge theories (LGTs) at non-zero temperature have been the focus of many investigations 
in recent years. Their study enables us to make non-perturbative predictions for the high temperature 
properties of gauge theories. Some of the issues which LGTs are expected to clarify are: the nature of the 
high temperature phase, the order of the phase transition, and the relevant elementary excitations at high 
temperatures. The pioneering work in was the first non-perturbative demonstration that quarks are 
dcconfmcd at high temperatures. This calculation is done in the strong coupling limit of the SU (2) LGT. 
Early Monte-Carlo simulations Q provided further support for the existence of this phase transition. Since 
then, there have been many studies of the thermodynamic properties of the SU(2) and the SU(3) LGTs 
H which have provided valuable insights into the high temperature behaviour of gauge theories. 

In this paper, we study the properties of the U(l) LGT at non-zero temperature. More precisely, we study 
the f7 (1) LGT on asymmetric lattices for which N a >> N T . N a is the lattice size in the spatial direction and 
N T is the lattice size in the temporal direction. The limit N„ — > oo and N T fixed, with periodic boundary 
conditions in the temporal direction, is the thermodynamic limit for a finite temperature system. There 
are several reasons why we have embarked on a study of this simple model. Firstly, the zero temperature 
properties of this model have been inferred from simulations on large symmetric lattices (N a = N T ) in 
In these simulations, a phase transition is observed from a confining phase to a deconfining phase. Though 
there was some controversy about the order of the phase transition, recent simulations on large lattices 
strongly suggest that the bulk transition is of first order j^]. The mechanism of this bulk transition is 
quite well understood in terms of monopole excitations [jll). These monopoles are topological objects that 
arise because of the periodicity properties of the action. The bulk system exists in two phases, a confining 
phase in which the monopole currents condense causing complete Meissner effect, and a deconfining phase 
in which the monopoles are too heavy to have any physical effect. It would be interesting to see how this 
picture of confinement vs deconfinement gets affected at finite temperature. This is also of relevance from 
the point of view of some recent claims of there being a second order phase transition in an extended 
version of the U(l) LGT with a monopole chemical potential term |lq| . A definition of a continuum theory 
is possible at the point of the second order phase transition, and an investigation of its properties can be 
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quite instructive. 

Secondly, the critical behaviour of gauge theories at high temperatures can be understood from analogous 
behaviour of three-dimensional spin models. The strong coupling analysis in jlj shows that the partition 
function of the SU (2) LGT can be rewritten as a three-dimensional spin model with a global Z{2) symmetry. 
The deconfining phase corresponds to the symmetry breaking phase of the spin model, and the confining 
phase corresponds to the symmetric phase. This symmetry arises because finite temperature gauge systems 
have an additional global symmetry arising from the periodic boundary conditions in the temporal direction. 
For systems with a gauge symmetry group G, this global symmetry group consists of elements belonging to 
the center of the gauge group. Furthermore, as emphasized in Q], the order of the phase transition in four 
dimensions is expected be dictated by the universality classes present in 3-d spin models having this global 
symmetry. These expectations have been borne out for the SU(2) || and the SU(3) LGTs |J in which one 
observes a second order Ising like and a first order Z(3) like phase transition respectively. One of our aims 
is to examine the same issue in the much simpler U(l) LGT. Unlike the non-abelian SU(N) LGTs, which 
have a discrete center subgroup, (Z(N)), the abelian U(l) LGT has a continuous center subgroup which 
is identical to the group itself. The role of the U(l) group on the confinement to deconfinement transition 
is also of some interest from the point of view of the abelian dominance hypothesis for confinement which 
holds that a U(l) subgroup controls the non-perturbative dynamics of non-abelian gauge theories |l9| ]. 
Finally, as the U(l) LGT already has a bulk phase transition unlike the SU(2) and SU(3) LGTs, there 
is the question of the interplay between this transition with the expected finite temperature transition. 
This matter has been recently examined in the context of SU(2) LGTs using a mixed action; it was found 
in ||lof that the bulk transition and the finite temperature transition may coincide, making it difficult to 
distinguish one from the other. Since simulations are always done on finite lattices, this raises the question 
whether the earlier studies on symmetric lattices pij were seeing a bulk transition or a finite temperature 
transition. As argued in p0[ , since "zero temperature" properties are inferred by extrapolating the results 
on symmetric lattices, and "finite temperature" properties are inferred by extrapolating the results on 
asymmetric lattices, it is not very easy to decide whether one is seeing a bulk effect or a genuine finite 
temperature effect. A major difference between a bulk transition and a finite temperature transition, which 
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can in principle distinguish the two transitions, is the movement of the transition point with the temporal 
lattice size. For non-abelian gauge theories, a physically relevant finite temperature transition should move 
towards the weak coupling region in a manner which is specified by the beta function of an asymptotically 
free gauge theory. With these notions in mind, we would like to investigate the finite temperature properties 
of abelian lattice gauge theories. 

We will mainly consider the Wilson action for the U(l) LGT, which is given by 

S = (3J2 cos(9(n pa/)). (1) 

The 9(n [iv) are the usual oriented plaquette variables: 

9(n iiv) = 6(n fj.) + 9(n + fi v) - 9(n + v fx) - 9(n v). (2) 

The link variables 9(n fj.) can take values from — ir to n. As mentioned before, the properties of this model 
at zero temperature are well known. There is a transition at (3 « 1.0 that is caused by a decondensation 
of monopole currents ( the bulk transition on a 16 4 lattices has been located precisely in pjJ and occurs 
at f3 = 1.016). These monopole currents are defined on the dual lattice by counting the number of Dirac 
strings entering or leaving a three dimensional cube on the original lattice pjj . The monopole density on 
a link (-kl) of the dual lattice is defined as 



<°N) = ^E »; (3) 

9{p) is extracted from 9(p) by expressing it as 

9(p) = 9(p) + 2irn(p), (4) 

so that it takes values from — tt to tt. In the above expression, the monopole current is defined on the links 
■kl of the dual lattice which are dual to the cubes c in the original lattice. As the monopoles form closed 
loops on the dual lattice, it is more convenient to measure the total perimeter density of the monopole 
current loops. The perimeter density is given by 
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Ni is the total number of links in the lattice. One observable which is of special relevance to our analysis 
is the Wilson-Polyakov line (henceforth called the Wilson line), which is defined as 

L (ft) = II exp(i0(n + n o 4 4)). (6) 

"0=0 

This observable is gauge invariant and has played a crucial role in studies of the deconfinement transition 
in SU (N) LGTs. The center symmetry that emerges at finite temperature is the transformation which 
multiplies by a constant phase all the time-like links emanating from some fixed time slice, namely 

exp(i8(n 4)) — ► exp(i a) exp(z 6{n 4)). (7) 

Although the action is invariant under this transformation, the Wilson line transforms as: 

L(n) — ► exp(i a)L(n). (8) 

Hence, a non-zero expectation value of the Wilson line signals a spontaneous breakdown of the center 
symmetry. The Wilson line is given the usual physical interpretation by writing it in the form: 

(L(n)) = exp(-/3F,(n)). (9) 

It measures the free energy (F q (n)) of a static charge in a heat bath at a temperature Hence, a non- 
zero value of (L(n)} indicates deconfinement of static charges whereas a zero value indicates confinement. 
We point out a subtlety in this representation, which has also been noticed before in the context of the 
SU(2) LGT The left hand side is in general a complex quantity and hence cannot be given a free- 

energy interpretation. One way of avoiding this problem is to work with the correlation function of two 
Wilson lines which is always a non-negative quantity. Apart from the Wilson line, we have also studied 
the susceptibility of the Wilson line which is defined as 

X i=N^)-{\g\n . (io) 

The " spin" variable s is constructed in Eq. [li]. The peak in the Wilson line susceptibility can be used to 
locate the phase transition. 

A related observable that is also relevant is the plaquette susceptibility, defined as 
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Xp = 67VXE« p2 >-^> 2 ) • (11) 
v 

In the above equation, P is the average plaquette density in the system. The peak in the plaquette 
susceptibility can also be used to locate the phase transition. According to the finite size scaling theory, 
the susceptibility of an observable, which blows up at a phase transition transition point in the infinite 
volume limit, is expected to scale with the system size as 

X c p r = N a ; (12) 

here N is the size of the system. For a first order transition, a equals d, the dimensionality of the system 
; for a second order transition, a equals 7/z/, where 7 and v are the exponents for the susceptibility and 
the correlation length near the transition, which are given by 

X = (T-T c )-t i={T-T c )~ v . (13) 

We will use the above observables to study the phase transition in the abelian lattice gauge theory. 

This paper is organized as follows. Sec 2 contains the results of our numerical investigations of this 
model. In Sec 3 we make some comments on the mixed action U(l) LGT and compare and contrast it 
with that of the of mixed action SU (2) LGT. In Sec 4 we summarize our conclusions. 

II. THE 17(1) LGT AT NON-ZERO TEMPERATURE. 

In this section we present our numerical analysis of the model. The system is mimicked at finite tem- 
perature by working on an asymmetric lattice (N a » N T ) with periodic boundary conditions in the 
Euclidean time direction. We first briefly describe the numerical procedure that we adopted to obtain our 
results. The Metropolis algorithm was used to generate successive Monte-Carlo configurations. A new link 
variable 9' was generated from the old one 9 by adding a number randomly chosen in the range (—a, a) 
with uniform probability. The value of a was tuned to get an acceptance of 50 percent. Care was taken 
so that the link variables remained in the range (— tt,tt). Simulations were performed on temporal lattices 
with N T = 2, 3 and 4. The observables that were measured are: the monopole density, the Wilson line, 
the plaquette density, the Wilson line susceptibility, and the plaquette susceptibility. As the Wilson line is 
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complex ( it is a phase with modulus equal to one ), we measure its real and imaginary parts separately. 
If we simply measure the average value of the real or the imaginary parts, the result will always be equal 
to zero because phase transitions are impossible on finite system as the tunnelling between the degenerate 
states always restores the symmetry. A rigorous way of studying symmetry breaking is to study the average 
value of the Wilson line in the presence of a small symmetry breaking external field, and then take the 
limit of zero field after taking the large volume limit. A simpler prescription, that is often employed in 
studying continuous spin models, is to study the root mean square value of the order parameter, and this 
is how we will proceed. The observable that we have measured as an indicator of spontaneous symmetry 
breaking is y/(s); s is defined as 



Here ReL and ImL refer to the average of the real and imaginary parts of the Wilson line respectively. A 
simple strong coupling analysis (valid for [3 small) yields the following effective action for the Wilson lines: 



N T is the temporal extent of the lattice and the 0(n) variables are the sums of the phases of all the time 
like links at the spatial point n. This is the action for the three dimensional planar model which is known 
to have an order-disorder transition at (3 cr — 0.454. For an N T = 2 lattice, this gives the critical coupling 
to be approximately 0.95. Thus we expect our lattice model to have a phase transition at f3 w 0.9. This 
strong coupling argument is valid only if the phase transition takes place in the strong coupling regime. 
Nevertheless, the strong coupling approximation provides a simple way of seeing how a three-dimensional 
spin model emerges from the four dimensional gauge theory 

A more direct way of seeing the appearance of an effective 3-d spin model, without using a strong 
coupling approximation, is to use the dual representation of the 4-d 17(1) LGT. The dual representation is 
given by 



This describes a gas of closed monopole loops (m M (r)) which interact by a 4-dim Coulomb potential. This 
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s = ReL 2 + ImL 2 . 



(14) 




(15) 




(16) 



system has a phase transition which takes place as a result of the competition between the energy and the 
entropy of the monopole loops. The fields </y(r) can be integrated out to give 



At non-zero temperature, the time direction is finite, and the Green's function satisfy periodic boundary 
conditions. The four dimensional Green's function can be rewritten as: 



here G M u (f — f") is the 3-dimensional Green's function. As the gas is finite in one direction, one again 
has a gas of monopole loops which now effectively interact with a 3-dim Coulomb interaction. From the 
point of view of the effective three-dimensional planar model, the monopole loops behave like vortex lines. 
The entropy of large loops in three-dimensions is smaller than in the four-dimensional case and this shifts 
the transition. The order of the transition cannot be deduced from these energy arguments and has to 
be determined by doing a finite size scaling analysis. We show in Fig. [I] the variation of y/{s) with j3 on 
a 6 3 2 lattice. The observable \J~{s(n)) is close to zero at small (3 and rises smoothly across the critical 
value. The t/(l) monopole density variation is shown in Fig. ^. There is a fall in the monopole density 
across the transition which coincides with the rise in the order parameter. In both cases, the variation is 
quite gradual and we would suspect that we are in the vicinity of a second order transition. This is made 
further suggestive by the gradual rise in the plaquette expectation value (see Fig. ^|). To determine the 
order of the transition, we perform a finite size scaling analysis of the susceptibility of the order parameter 
(xz)- We have done the finite size scaling study on lattices of temporal extent N T = 3 and N T — 4. For 
this method to work, it is crucial that we are very close to the pseudo-critical point corresponding to the 
lattice that we are working on. The histogram method is used to extrapolate observables from one value 
to a nearby neighbouring value. The pseudo-critical point is located in this way by looking at the peak in 
the susceptibility of the order parameter. The behaviour of the susceptibility near the transition on 6, 9, 12 
and 16 sized spatial lattices (keeping the temporal extent fixed at N T — 3) is shown in Fig. ^. These results 
were got after 200,000 measurements on N a = 6,9 lattices while 150,000 measurements were made on the 
N a — 12 lattice and 100,000 measurements were made on the N a — 16 lattice. All the measurements were 
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(17) 



(18) 
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made after ignoring the first 50,000 Monte-Carlo iterations. The errors were estimated by binning the 
data. The finite size scaling theory predicts 

X ~ Nt (19) 

for a first order phase transition; d being the effective dimension of the system (in this case d — 3); for a 
second order transition, the prediction is 

X~Nj; (20) 

7 and v are the exponents in Eq. [l3|. The peaks in the susceptiblity were fit to an N" dependence and 
a good fit is obtained (Fig. ||). From our fit we determine <y/v to be 1.95 with an error of 0.13. For the 
three dimensional planar model, which is the effective spin model with which we would like to compare our 
results, the ratio j/u = 1.97. 

We have performed another test for the order of the phase transition by studying the distribution of 
the monopole loops. Since the monopole loops, which now behave as vortex lines, are responsible for 
driving the phase transition, they should also exhibit a scaling behaviour at the point where the transition 
becomes second order. At the point of the second order phase transition, there are fluctuations on all 
length scales and the monopole loops come in all sizes and shapes. At the transition point, the monopole 
loop distribution function, p(l), should scale as 

P(l) = p (21) 

with an exponent r which is independent of the volume. The function p{l) is defined as the probability 
of finding a loop of length I at a site. We have calculated p(l) for various lattices at their pseudo-critical 
points and we find that p(l) can be nicely fitted to the above form. The value of r that we get (using 
data on a N a — 16 lattice) is 2.27 ± .003. Unfortunately, we are not able to compare it with a theoretical 
calculation. Nevertheless, it is interesting that the topological excitations, in this case the monopoles, 
exhibit a scaling dependence at the transition point. Fig. ^ shows p(l) for loops of length I upto 14 on 
different lattices. The plot shows that they fall on a straight line independent of the volume. The slight 
distortions for the N a — 6 lattice are due to finite size corrections when the loop length is larger than the 
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lattice size. A study of larger loop lengths can also be made but the algorithm that we have used to count 
the loops slows down when the loop lengths become very large, and so we have not studied larger loops. 
Away from the critical point, p(l) either falls rapidly to zero (in the deconfining phase) or remains constant 
(in the confining phase). A similar study of the monopole loop distribution function was made in p2j for 
the abelian projected monopoles in the SU(2) theory near the continuum limit. 

Now we say a few words on the location of this transition as compared to the location of the bulk 
transition. The transition point (which is located by looking at the peak in the susceptibility of the order 
parameter) shifts as a function of N T as: 

N T 2 3 4 

P cr 0.9297(3) 1.012(2) 1.032(2) 

The above critical values are those obtained on N a of 16, 16 and 12 respectively. The bulk transition on a 
16 4 lattice was located at /3 = 1.016 in f§. Since the bulk transition is known to be of first order and 
the transition that we have observed is of second order, we are observing a change in the order of the phase 
transition as a function of N T . The nature of the transition that we have observed can be further studied 



by monitoring the behaviour of the plaquette susceptibility which was defined in Eq. 11. The behaviour of 
the plaquette susceptibility for a 6 3 3 and a 12 3 3 lattice is shown in Fig. [| This graph clearly shows that 
the plaquette susceptibility does not scale with the four dimensional volume as is expected of a first order 
bulk transition. 

The above observations show that the transition that we have observed (on an N T = 3 lattice) is a 
dcconfincmcnt transition whose scaling behaviour is quite distinct from that of the bulk transition. The 
order of the transition as seen by the Wilson line is a second order transition with the ratio j/v which 
has approximately the same value as that in the three-dimensional planar model. Since the plaquette 
susceptibility does not scale as a first order bulk transition, we are not observing the bulk transition at 
the point where we are observing the finite temperature transition. The finite size scaling analysis can be 
repeated on an N T — 4 lattice and the peak in the susceptibility of the order parameter as a function of 
volume is shown in Fig. |^. This scaling again suggests a second order transition with a value for j/v which 
is again in good agreement with that of the three-dimensional planar model. 
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We now place our results in the context of some recent developments. The studies in [[L0| considered the 
finite temperature properties of the mixed action SU (2) LGT which is defined by 

S = {fi f /2) tr f U(p) + 0S„/3) £ tr a U(p). (22) 
p v 

This model is known to have lines of first order bulk transitions in the (3f,f3 a plane p3]. It was found in 
fiof that the deconfinement transition of the pure ST/ (2) LGT, which is of second order, continues into the 
Pffia plane and joins the line of first order bulk transitions. Based on this observation in pj, a possibility 
was considered where the there is either only a bulk or a finite temperature transition in the SU (2) LGT. 
On the other hand, there may also be a very small but finite separation in the two transitions which cannot 
be resolved on the lattice sizes used in the simulations. In our case, the relevant coupling space is the j3, T 
plane. Our simulations of the U(l) LGT show that there is a deconfinement transition which is of second 
order and that the plaquette susceptibility does not scale as a first order bulk transition at the transition 
point. The transition is shown to shift from its bulk value as a function of N T but this shift is very small, 
though still discernible. We also place our results in the light of the expectations of Svetitsky and Yaffe 
. According to their general arguments, the critical begaviour of the £7(1) LGT at finite temperature is 
expected to fall in the same universality class as that of the 3-d planar model, which is known to have a 
second order phase transition. Our finite size scaling analysis on the N T = 3 lattice definitely rules out a 
first order phase transition and indicates that the ratio 7/1/ has approximately the same value as in the 
3-d planar model. 

We conclude this section with a proposal for the finite temperature phase diagram for the U(l) LGT. 
In the bulk system (which corresponds to taking the limit N„ = N T and N a — > 00), there is the monopole 
driven phase transition. On asymmetric lattices of very small temporal extent, we expect the bulk system to 
look like a three dimensional £7(1) gauge theory. The three-dimensional £7(1) LGT has no bulk transition. 
On asymmetric lattices of very large temporal extent we expect to see the bulk transition. Lattices of 
intermediate temporal extent will exhibit a complicated crossover from a four-dimensional bulk system to 
a three-dimensional bulk system. Our simulations on an N T = 3 lattice show that there is a second order 
phase transition which has the same value of •y/v as in the three-dimensional planar model. This transition 
is a deconfinement transition as is indicated by the behaviour of the fundamental Wilson line. However, 
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there is no bulk transition on this lattice. 



III. MIXED ACTION 17(1) LGT. 

Since lattice actions are anyway not unique, we can always construct more complicated looking lattice 
actions and examine their properties. A simple generalization of the action in Sec 1 is the mixed action 
which is defined by 

S = ^J2 cos(0(p)) + 2 cos(20(p)). (23) 

The two pieces of the above action are different from each other only insofar as their periodicity properties 
are concerned. In the naive continuum limit, a — ► 0, the second term is like an irrelevant coupling and 
is not expected to change the long distance properties of the theory. The zero temperature properties of 
this action have been studied in Jl7t and it has a rich phase structure of first and second order transitions. 
This system can also be studied at finite temperature just as we did for the 0% = theory. Again, from a 
simple strong coupling analysis of the mixed action £7(1) LGT, we get an effective theory of spins which is 
that of the mixed planar model 

S eff = 2(§0^ X> S (0(n) - 0{n')) + 2(|)^£ cos(20(n) - 26{n>)). (24) 

nn' fin' 

Putting Pi = gives a three-dimensional planar model of spins, the only difference being in the periodicity 
properties of the action. Hence, the finite temperature properties of the mixed model in the Pi = limit 
should be identical to those in the P% — limit. In particular, the previous statements regarding the 
order of the transition will also be true in this limit. A surprising feature of the mixed planar model is 
that it posesses a region of first order phase transitions for some values of P2 p6[ | . This implies a similar 
region of first order transitions in the mixed U(l) LGT for a segment of P2 values. The order of the finite 
temperature transition changing in the direction of an irrelevant coupling has also been discussed in the 
context of the mixed action SU(2) LGT 0. 

From the above analysis, it is clear that there are many similarities between the C7(l) LGT and the 
mixed action SU(2) LGT. We now show that the similarity also extends to the construction of the relevant 
order parameters in the two models. The mixed action SU(2) LGT Q is defined by 
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s = y E WW + y E ^(p); ( 25 ) 
p p 

iry and ir a denote the traces in the fundamental and the adjoint representations respectively. The limit 
Pa — describes an SU(2) LGT and the limit (3f = describes an SO(3) LGT. The order parameter of 
the finite temperature transition in the SU (2) LGT is the Wilson line in the fundamental representation 
which is defined as 

L f (n) = Tr f ]J U{n + n 4 4). (26) 

710=0 

In the 5*0(3) LGT (j3f — 0), this observable is identically zero because of the following local Z(2) symmetry: 

U(n 4) -> Z{n)U{n 4); (27) 

Z(n) can take the values +1 or —1 at any site. For the SO(3) LGT, the appropriate order parameter is 
the Wilson line in the adjoint representation, 

L a (n) =Tr a jj U(n + n a 4 4) , (28) 

"0=0 

which is invariant under the local Z(2) transformation in Eq. |2^. For the group SU(2), Lf and L a are 
related by 

L a (n) = L f (nf - 1.0 . (29) 

In the mixed action U(l) LGT we are faced with a similar problem in defining the order parameter for the 
Pi = theory. In this limit, the Wilson line defined in Sec 1 is identically zero because of the following 
local symmetry: 

exp(i6»(n 4)) -> Z(n) exp(i 6 (ft 4)). (30) 
The correct order parameter to use in this limit is 

L2(n) = n exp(i 20(n + n o 4 4)) , (31) 

Tl =0 

which is analogous to the Wilson line in the adjoint representation of SU(2). The relationship between 
Li2(n) and L(n) is Li2(n) = Lin) 2 . 
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Though there are similarities at the formal level, and in the phase structure of the mixed action abelian 
and mixed action non-abelian theories, there are, of-course, many important differences between the two 
systems. An important difference is that, unlike in the abelian LGT, the critical coupling in the SU(2) 
LGT is expected to scale with N T according to the beta function of the Yang-Mills theory. There is strong 
evidence for this asymptotic scaling from simulations on very large temporal lattices p5| ] . 

The purpose of this section was only to indicate that many of the issues such as the mixing of the 
bulk and finite temperature transitions which have been raised recently can all be explored in the much 
simpler mixed action abelian lattice gauge theory. Also, because the physical properties of this model are 
well understood in terms of the monopole excitations, this model may prove useful in investigating these 
issues. 

IV. CONCLUSIONS 

In this paper we studied the U(l) LGT using the Wilson line as the order parameter. We have found that 
there is a transition into a deconfining phase at large coupling which is driven by the decondensation of 
monopole loops. The monopole loops, however, effectively interact like the vortices in the three-dimensional 
planar model. The deconfining phase breaks the global U(l) symmetry present in the theory. A finite size 
scaling analysis of the susceptibility of the Wilson line indicates that the transition is of second order ( 
on lattices of temporal size N T = 3 and N T = 4) with a ratio for j/v which has almost the same value 
as in the three-dimensional planar model. A scaling form for the distribution of the monopole loops was 
also established at the point of the second order phase transition. This transition was also examined by 
studying the plaquette susceptibility at the transition point. The plaquette susceptibility (on a N T = 3 
lattice) does not scale as is expected of a first order bulk transition. There is also a small shift in the 
transition point from the bulk value. We have also pointed out that many of the recent issues concerning 
the mixing of the bulk and finite temperature transitions can also be raised in the abelian lattice gauge 
theory Since the abelian theories are well understood in terms of their monopole excitations, some of these 
issues can perhaps be clarified. 
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1G 



FIG. 1. Variation of the order parameter on a 6 3 2 lattice. 
FIG. 2. Plaquette expectation value on a 6 3 2 lattice. 
FIG. 3. Monopole density on a 6 3 2 lattice. 
FIG. 4. Susceptibility of the order parameter near the transition on 6, 9, 12 and 16 size spatial lattices. 

FIG. 5. A fit of the maximum value of the logarithm of the susceptibility to InN. 

FIG. 6. Scaling of the monopole loop distribution function p(l) with loop length (1). The graph shows the 
logarithm of p(l) as a function of the logarithm of I. 

FIG. 7. A finite size study of the susceptibility of the order parameter on a N T = 4 lattice. 

FIG. 8. The scaling of the plaquette susceptibility near the phase transition. This is shown on 6 3 3 and 12 3 3 
lattices. 
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